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Control of atomic currents using a quantum stirring device 
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PACS 03 . 65 . -w - Quantum mechanics. 

PACS 03.65.Vf - Phases: geometric; dynamic or topological. 
PACS 05 . 30 . Jp - Boson systems. 

Abstract. - We propose a BEC stirring device which can be regarded as the incorporation of a 
quantum pump into a closed circuit: it produces a DC circulating current in response to a cyclic 
adiabatic change of two control parameters of an optical trap. We demonstrate the feasibility 
of this concept and point out that such device can be utilized in order to probe the interatomic 
interactions. 
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The realization of Bose-Einstein Condensation (BEC) 
of ultra-cold atoms in optical lattices and atom chips [1] 
and the availability of conveyor belts [2] is considered to be 
a major breakthrough with potential applications in the 
arena of quantum information processing [3], atom inter- 
ferometry [4-6], lasers [1,7] and atom diodes and transis- 
tors [8-10]. A major advantage of BEC based devices, as 
compared to conventional solid-state structures, lies in the 
extraordinary degree of precision and control that is avail- 
able, regarding not only the confining potential, but also 
the strength of the interatomic interactions, their prepa- 
ration and the measurement of the atomic cloud. Accord- 
ingly it is envisioned that the emerging field of "atomtron- 
ics" will provide a new generation of nanoscale devices. 

The possibility to induce DC currents by periodic (AC) 
modulation of the potential is familiar from the context 
of electronic devices. If an open geometry is concerned, it 
is known as "quantum pumping" [11] , while for closed 
geometries we use the term "quantum stirring" [12]. We 
consider below the stirring of condensed ultra-cold atoms 
due to the periodic variation of the on-site potentials 
and of the tunneling rates between adjunct confining 
traps. We show that the nature of the transport process 
depends crucially on the sign and on the strength of 
the interatomic interactions. We distinguish between 
four regimes of dynamical behavior: For strong repulsive 
interaction the particles are transported one-by-one, 
which we call sequential crossing; for weaker repulsive 
interaction we observe either gradual crossing or coherent 





Fig. 1: (Color online) Illustration of the model system. In 

the first half of the cycle a the particles are transported from 
the gated-well to the double-well via the fci bond, while in the 
second half of the cycle b the particles are transported back 
from the double-well to the gated-well via the ki bond. See 
the text for further details. 



mega crossing; finally, for strong attractive interaction 
the particles move together as one composite unit from 
trap to trap. 

Model — The simplest model that captures the physics 
of quantum stirring is the three-site Bose-Hubbard Hamil- 
tonian (BHH) [10, 13-15] (see Fig. 1) 1 . We call the i=0 
site "gated-well" since we late assume that we have con- 
trol over its potential energy vq = e. The i—1,2 sites form 



1 The three-site BHH is a prototype system for many recent stud- 
ies [13]. A classical analysis of this system was performed in [15] 
where it was shown that for appropriate system parameters and ini- 
tial conditions chaotic dynamics would emerge. In this work we 
consider adiabatic driving of ground state preparation and therefore 
chaotic motion is not an issue. 
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a two level "double- well" with potential energy vi=v 2 = 
The N boson BHH is: 



2 r 



v i n l +—n l (n i -l) 



-kMpi- 



^fc i &jSo+h.c. (1) 



We set h=l which corresponds to measuring energies in 
units of frequency. Furthermore, without loss of general- 
ity we choose time units such that k c —l. Accordingly the 
two single particle levels of the double- well are e± — ±1. 
The annihilation and creation operators Sj and b\ obey 
the canonical commutation relations [Si, 6^] = Si.j while 

the operators hi = b\bi count the number of bosons at 
site i. The interaction strength between two atoms in a 
single site is given by U — Aita s V B{{ /m where V off is the ef- 
fective volume, m is the atomic mass, and a s is the s-wave 
scattering length which can be changed by applying an 
additional magnetic field [16]. The on-site potentials Vi as 
well as the coupling strengths fcj are controlled by changing 
the confining potential. The main assumption underlying 
the BHH is that the single-particle ground state wave- 
functions are sufficiently localized at the sites, and that 
for the temperatures involved they are well separated in 
energy from the excited single-particle levels. Experimen- 
tally, such deep trapping potentials can accommodate sev- 
eral hundred particles [17]. 

The couplings between the gated-well and the two 
ends of the double- well are k\ and k 2 . We assume that 
both are much smaller than k c (for the two-mode BEC 
dynamics see for example [18]). It is convenient to 
define the two control parameters of the pumping as 
X\ = (l/fc 2 ) — (l/&i) and X 2 — e. By periodic cycling of 
the parameters (X\, X 2 ) we can obtain a non-zero amount 
(Q) of transported atoms per cycle. The pumping cycle 
is illustrated in Figs. 1,2. Initially all the particles are 
located in the gated-well which has a sufficiently negative 
on-site potential energy (X 2 < 0). In the first half of 
the cycle the coupling is biased in favor of the k\ route 
(Xi > 0) while X 2 is raised until (say) the gated-well is 
empty 2 . In the second half of the cycle the coupling 
is biased in favor of the k 2 route (Xi < 0), while X 2 
is lowered until the gated-well is full. Assuming U=0, 
the gated-well is depopulated via the k\ route into the 
lower energy level e_ during the first half of the cycle, 
and re-populated via the k 2 route during the second 
half of the cycle. Accordingly the net effect is to have 
a non-zero Q. If we had a single particle in the system, 
the net effect would be to pump roughly one particle 
per cycle. If we have N non-interacting particles, the 
result of the same cycle is to pump roughly N particles 
per cycle 3 . We would like to know what is the actual 



2 Later we estimate the required X2 variation in order to have all 
the particles transferred from the gated-well to the double well. We 
also analyze smaller pumping cycles (see Fig. 3b) for which only a 
fraction of particles gets through. 

3 Through the driving cycle the total number of bosons remains 
constant. The energy is not a constant of motion, but in the adia- 
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Fig. 2: Illustration of the pumping cycle. See the text for fur- 
ther details. For a large cycle that encircles the whole shaded 
region we have Q w N. Position of the monopoles: a no in- 
teractions, b with interactions, c, the energy levels along the 
Xi = axis are schematically plotted for an N = 3 system. 



result using a proper quantum mechanical calculation, 
and furthermore we would like to investigate what is the 
effect of the interatomic interaction U on the result. 

Methods — Within the framework of linear response 
theory the induced current is / = — GxXi if we change Xi 
and I = —G 2 X 2 if we change X 2 . The coefficients Gi and 
G 2 in these linear relations are defined as the elements of 
the geometric conductance matrix, and can be calculated 
using the Kubo formula approach (see below) . Integrating 
the current over a full cycle we get 



Q= <j> Idt= - <j)(G 1 dX 1 + G 2 dX 2 ). 



(2) 



In order to calculate the geometric conductance we use the 
Kubo formula approach to quantum pumping [19] which 
is based on the theory of adiabatic processes [20] . It turns 
out that in the strict adiabatic limit G is related to the 
vector field B in the theory of the Berry phase which is 



batic limit the system comes back to the same state at the end of 
each cycle. 
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known as the "Berry Curvature" . The adiabatic slowness 
condition on X in the present context, taking into account 
the two-orbital approximation, is discussed in Section 4 
of [21]. Using the notations £?i = — Gi and B 2 = G\ it is 
illuminating to rewrite Eq.(2) as 



= <h B ■ ds, 



(3) 



where we define the normal vector ds — (dA^,— dX{) as 
illustrated in Fig. 2. The calculation of the so-called Kubo- 
Berry Curvature is done using the following formula [19]: 



B: 



E 



2 Im[Z~ non ] J~ n7 



Above I = i/2 fci(6 6i-&|fe ) + k 2 (blb -blb 2 ) 
eraged current 4 via the bonds Oi — >1 and 2 



E ) 2 



(4) 

is the av- 
-+0, while 

= —dH/dXj is the generalized force associated with 
the control parameter Xj. The index n labels the eigen- 
states of the many-body Hamiltonian. We assume from 
now on that no is the BEC ground state 1 . 

The advantage of the above, so called "geometric" 
point of view is in the intuition that it gives for the 
result: Formally the field B is like a projection of a fic- 
titious magnetic field in an embedding three-dimensional 
X space (the third coordinate A3 if formally defined 
such that X = —dH/dX$). The flux of this fictitious 
magnetic field through any out-of-plane surface which is 
enclosed by the pumping cycle gives the so-called Berry 
phase, while the line integral over this fictitious magnetic 
field, i.e. Eq.(3) 5 , gives Q. As implied by inspection of 
Eq.(4) the sources of B are located at points where the 
ground level no has a degeneracy with the next level. 
A simple argument implies that this "magnetic charge" 
is quantized like Dirac monopoles, else the Berry phase 
would be ill-defined. For details see [19]. In our model 
system for U—0 all the "magnetic charge" is concentrated 
in one point. As the interaction U becomes larger 
the (AT+l)-fold degeneracy of the levels is lifted, and 
this "magnetic charge" disintegrates into N elementary 
"monopoles" (see Fig. 2). We further discuss the energy 
spectrum in the next section. 

Regimes — We define the average coupling as n — 
(ki+k 2 )/\/2. In the zeroth-order approximation k\ and 
k 2 are neglected, and later we take them into account as 
a perturbation. For k=0 the number (n) of particles in 
the gated-well becomes a good quantum number hence we 
can associate the level index n with the number of particles 

4 Due to the continuity equation the same result for Q is obtained 
irrespective of whether we measure the current via the Oi — »1 bond or 
via the 2i— >0 bond. The advantage of using the the "symmetrized" 
version for X is that the same amount of particles is being transported 
during both halves of the cycle, allowing us to focus on (say) the first 
half and then double the result. 

5 In the X$=0 plane the third component of the fictitious magnetic 
field is zero due to the time-reversal symmetry [19]. 




Fig. 3: (Color online) Conductance during the first half of 

the pumping cycle, a, the conductance Gi as a function of 
the on-site potential e, for various values of U. In the nu- 
merical simulations (see also Fig. 4) we are using exact di- 
agonalization of the trimer Hamiltonian (1). For the evalu- 
ation of G we use Eq. (2) while for Q we use Eq. (4). The 
other parameters are N = 16 particles, « = 0.0003/^ and 
(ki — fe)/v2 = 0.0001/\/2. As the interaction U becomes 
larger one observes the crossover from a single to individual 
peaks in the conductance, b, the [/-dependence of the inte- 
grated charge Q* , calculated for wide rectangular cycles for 
which X2 is varied within [— 00, £*]. The values of e* are indi- 
cated by arrows of the same color in the main panel. 



in the gated-well. Furthermore, we adopt a "two-orbital 
approximation" : we assume that there is non-zero occu- 
pation only in the gated-well and in the lower double- well 
level, which is valid if NU -C k c . Note also that we assume 
an adiabatic process 3 , and accordingly non-adiabatic tran- 
sitions to the higher orbitals can be safely neglected [21]. 

Within the two-orbital approximation the many body 
energies are E„ = £' gatcd _ wc „(n) + E doublo _„ cll (N~n), where 
n=0,l, -,N, and E g&ted _„ eU (n) = en + (l/2)U(n— l)n, and 
E dmblQ _ mll (N~n) = -(N-n) + (l/4)U(N-n-l)(N-n) . 
The location e n = -1 + \U x (N-in+2) of the 
n 1 ► (n— 1) crossing is determined from the 
degeneracy condition E n — B n _i = 0, where 
n=l,2, ...N. The N crossings are distributed within 
e € -1 + (N-l)[-U, U/2]. The rescaled version of the 
control variable is e = (e+l)/((N— l)U), and its support 
is — 1 < e < 1/2. The distance between the crossings, 
while varying the gated-well potential e, is U. Once we 
take k into account we get avoided crossings of width 
5e n = [(N+l— n)n/2] 1 / 2 n. If k/UN is large these avoided 
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crossings merge and eventually we get one mega crossing. 

For the purpose of further analysis we apply the 
two-orbital approximation, within which the many 
body Hamiltonian matrix is H nm — E n 5 n . m — K n 5 n ^ n ±\, 
where n=Q, ■•• ,N and the couplings are defined as 
n n = (n— l\7i\n). The calculation involves the matrix 
elements of b\bo, leading to K n = [(N + 1 — ra)™] 1 / 2 n. 
Analogous expression applies to the current operator 
where k is replaced by (fci — k 2 )/\/2. For large U, as e is 
varied, we encounter (say for N = 3) a sequence of distinct 
Landau-Zener transitions (|3) i— > |2) |1) |0)). The 
distance between avoided crossings is of order U while 
their width is Ss n — K n . The widest crossings are at the 
center with Se n <~ Nk. This should be contrasted with 
the energy scales U and NU that describe the span of 
the crossings. Accordingly we deduce that for repulsive 
interaction there are three distinct regimes: for U <C k/N 
we have a "mega crossing" ; for U 3> Nk we have the 
"sequential crossing regime", while in the intermediate 
regime we have a "gradual crossing" . Below we summa- 
rize the results in the various regimes. In particular Eq.(7) 
is obtained from Eq. (4) with a two level approximation for 
each crossing. We also related briefly to the U < regime. 

Results — Most of the contribution to the line integral 
in Eq.(2) comes as we change X 2 during the avoided cross- 
ings that have been discussed in the previous section. If we 
close the pumping cycle outside of this limited range, then 
the Xi-variation can be safely neglected. Accordingly we 
refer from now on to G 2 = G only. An overview of the 
numerical results for the conductance is shown in Fig. 3, 
where we plot G as a function of X 2 for various interaction 
strengths U. In the same figure we report the normalized 
amount of particles Q for various driving cycles. As the 
shape of G 2 changes, the dependence of Q on the X 2 -span 
of the pumping cycle becomes of importance. Thus, by 
measuring Q we obtain information on the strength of the 
interatomic interactions. 

In Fig. 4 more details are presented: besides G we also 
plot the A 2 -dependence of the energy levels, and of the site 
population. Four representative values of U are considered 
including also the U < case. Let us discuss the observed 
results. For U—0 all the particles cross "together" from 
the gated-well orbital to the e_ double-well orbital. We 
call this type of dynamics "mega crossing" . The outcome 
is just N times the single particle result: 



G 



-N- 



(5) 



(kl-kj)/2 
[( £ -e_) 2 + 2(fc 1 +fc 2 ) 2 ]3/ 2 ' 

which can be expressed in terms of the control parame- 
ters (Xi,X 2 ). This result approximately hold as long as 
U <C k/N . Integrating over a full cycle one obtains 



Q = N 



[1 + ( K R) 2 } 1/2 - I 



kR 



(6) 



where R is the radius of the pumping cycle (see Fig. 2). 
For small cycles we get Q « NkR/2, while for large cycles 



we get the limiting value Q N. In the other extreme, 
for very repulsive interaction (U^Nk) we get 



G 



N 



fci - k 2 

h + k 2 J ^ [(e - e„) 2 + (2<5e„) 2 ]3/ 2 ' 



(7) 



For intermediate values U G [k/N, Nk], we find neither 
the sequential crossing of Eq.(7), nor the mega-crossing 
of Eq.(5), but rather a gradual crossing. Namely, in this 
regime, over a range AX 2 = (3/2)(iV— 1)U we get a con- 
stant geometric conductance: 



G 



fci 



fci + k 2 



1 

3C7 



(8) 



which reflects in a simple way the interaction strength. 
This formula was deduced by extrapolating Eq.7, and then 
was validated numerically (lower panel of Fig. 4c) . 

As discussed above for large positive U the (N + l)-fold 
"degeneracy" of the U =0 Landau-Zener crossing is lifted, 
and we get a sequence of N Landau-Zener crossings (for 
schematic illustration see the lower panel of Fig. 2, and 
compare with the numerical results in the upper panels 
of Fig. 4). Also for U<0 this (7V+l)-fold "degeneracy" is 
lifted, but in a different way: the levels separate in the 
"vertical" (energy) direction rather than "horizontally" 
(see upper panels of Fig. 4). Accordingly all the particles 
execute a single two-level transition from the gated-well 
to the double- well (see Fig. 4a) . This direct Landau-Zener 
transition from the n=N level to the n—0 level is very 
sharp because it is mediated by an 7V-th order virtual 
transition via the intermediate n states. Accordingly, for 
sufficiently strong attractive interaction all the particles 
move together from the gated-well to one of the double- 
well sites. When the sign of X\ is reversed they are 
transported from one end of the double-well to the other 
end (not shown). This should be clearly distinguished 
from the (N + l)-fold degenerated transition to the lower 
double-well level which is observed in the U=0 case. 

Summary — The theoretical [10, 18,22] and experimen- 
tal [23] study of driven dynamics in single and double 
site systems is the state of the art. Study of three-site 
systems adds the exciting topological aspect: controlled 
atomic current can be induced using optical lattice tech- 
nology [24]. The actual measurement of induced neutral 
currents poses a challenge to experimentalists. In fact 
there is a variety of techniques that have been proposed 
for this purpose. For example one can exploit the Dopplcr 
effect in the perpendicular direction, which is known as the 
rotational frequency shift [25] . The analysis of the proto- 
type trimer system reveals the crucial importance of inter- 
actions. The interactions are not merely a perturbation: 
rather they determine the nature of the transport process. 
We expect the induced circulating atomic current to be 
extremely accurate, which would open the way to various 
applications, either as a new metrological standard, or as 
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Fig. 4: (Color online) Evolution of the energy levels, the site occupation and the conductance. Further details relating to the 
data of Fig. 3. We refer to four representative values of U, which are indicated on top of each set of panels. Upper panels: the 
lowest N+l energy levels E n which dominate the conductance G2 are plotted as a function of ~K% = e. The insets represent 
magnifications of the indicated areas. Middle panels: the site occupations no (blue A), n\ (black o),n2 (red □). Lower panels: 
the corresponding conductance G2 as a function of e. Numerical results are represented by solid black lines while the dotted 
red line corresponds to the analytical result (5) in b and to (7) in c, d. 



a component of a new type of quantum information or 
processing device. 
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